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Our main aim in this paper is to promote the coframe variational method as a unified approach 
to derive field equations for any given gravitational action containing the algebraic functions of the 
scalars constructed from the Riemann curvature tensor and its contractions. We are able to derive a 
master equation which expresses the variational derivatives of the generalized gravitational actions 
in terms of the variational derivatives of its constituent curvature scalars. Using the Lagrange 
multiplier method relative to an orthonormal coframe, we investigate the variational procedures for 
modified gravitational Lagrangian densities in spacetime dimensions n 3. We study well-known 
' gravitational actions such as those involving the Gauss-Bonnet and Ricci-squared, Kretchmann 

scalar, Weyl-squared terms and their algebraic generalizations similar to generic f{R) theories and 
^ ' the algebraic generalization of sixth order gravitational Lagrangians. We put forth a new model 

' involving the gravitational Chern-Simons term and also give three dimensional New massive gravity 

, equations in a new form in terms of the Cotton 2-form. 
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I. INTRODUCTION 



Although General Relativity theory is an experimentally well-tested and successful theory [1], the search for al- 
ternative or modified gravity theories emerged immediately after the formulation of General Relativity (GR). The 
motivation of introducing such theories varies considerably. For example, in some of the earlier theories, such as those 
introduced by Kretschmann [2], Weyl [3] and Eddington [4], Lagrangians containing higher order curvature scalars 
^ were considered in order to understand mathematical and physical machinery of Einstein's theory and to compare 
fSl • their predictions with different and more generalized theories. In another important route, in the same era, Kaluza 

■ and Klein showed that it is possible to unify gravity with electromagnetism from pure geometric Einstein-Hilbert 
^vq ; (EH) action but employing an extra compact dimension [5] . Although this theory is not a viable theory, its main idea 
y—{ ' of employing extra dimensions for the unification survive in modern theories such as String/M theory [6]. In these 

theories, terms involving higher order curvature scalars arise naturally in gravitational sector of their action [7]. It 
is also known that, in order to incorporate quantum effects to gravity, one might employ Wheeler-DeWitt effective 
] action formalism [8], which leads to higher order curvature corrections as well [9]. Hence, the studies considering 
' ^ ' unification as well as quantization of gravity require modifications of some sort in the gravitational Lagrangian. 

■ Another motivation for considering modified gravitational actions comes from the dynamics of the universe. The 
^ . results of data accumulated on cosmology and astrophysics, in recent years, lead to interesting observations, for 
- - ' example, that the universe is in the epoch of accelerating expansion driven by the dark energy and most of the matter 

content of the universe is dark [10]. Although the ongoing and future particle physics experiments may explain the 
physical content of dark matter and dark energy, the possibility that these observations suggests that GR should be 
modified is also an interesting and attractive idea. Hence, in recent years, the old idea of replacing R in EH action 
with /(i?) [11-15], has been studied extensively in this context. For more information as well as recent developments 
and their applications to the physically relevant models of f{R) theories, see one of the excellent reviews [16, 17, 19] 
and references therein. The problem of this approach is that, can these modified gravity theories explain those 
observational phenomena without bringing out more problems to the area, such as instabilities, and can be compatible 
with laboratory and observational astrophysical data. Actually, one should regard these attempts as toy models, which 
might help to understand the gravitational aspects of nature by modeling and studying geometric gravity theories. It is 
preferable to start investigating modified gravitational models by starting from action scalar for which the metric field 
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equations derived from an action via a well-defined variational procedure. There are distinct variational procedures 
that lead to distinct field equations. In so-called Palatini formalism, the metric tensor and connection are treated 
as independent variables of the theory and the relation between them is derived as a consequence. In contrast, in 
the metric formalism, the metric tensor is assumed to be the only field variable determining the spacetime geometry 
and the field equations for it follows from the variational derivatives of the gravitational action with respect to the 
components of the metric tensor. Except for the Einstein-Hilbert Lagrangian the metric field equations that follow 
from these two procedure, as we shall explicitly show below, are not the same. In addition, in the metric approach, 
the connection (Riemannian or Levi-Civita connection) is derived uniquely from the derivatives of the metric tensor, 
whereas in the Palatini approach, the connection may be non-Riemannian. This has important consequences on the 
motions of test particles since a connection defines parallel transport and therefore determine the geodesic structure 
of the spacetime. 

The organization of our paper can be outlined as follows. The paper can roughly be divided into three parts. The 
first part deals with the constrained variational procedure and Lagrange multiplier technique that will be used in 
the rest of the paper. The second part is related to the algebraic generalization of modified gravitational actions. In 
this part, the variational derivative of the well-known f{R) theory in both metric and Palatini approaches and the 
equivalence of these actions to the Brans-Dicke (BD) type scalar tensor theories are discussed. In the subsequent 
sections we calculate the variational derivative of the actions which consist of functions of the terms such as Gauss- 
Bonnet, Ricci-squared or Weyl-squared terms and their algebraic generalizations. Some of the results related to the 
algebraic generalization of the quadratic curvature terms are new. Various important inter-relationships between the 
Lagrangians in three and four dimensions are indicated. Finally, the variation of terms involving derivatives of the 
curvature scalars, which lead to sixth order theories, as well as their algebraic generalizations are discussed. The third 
part is the original part and it consist of last two sections. These sections can be considered as applications of the 
general formulae obtained in the previous parts. In the first of these sections, namely section V, a new model for the 
Chern-Simons gravity in four dimensions is derived whereas the last section is devoted to the study of field equations 
of three dimensional general massive gravity models and in particular Lagrange multipliers of general massive gravity 
theory. Our notation and conventions are briefly explained in the Appendix. 



Wc; find it convenient to start our discussion in the next section with the general definitions and the remarks about 
the variational procedure we shall use throughout the paper. The field equations of Einstein's General Theory of 
Relativity can be derived from Einstein-Hilbert (EH) action 



via a variational procedure, where the gravitational Lagrangian is given by Ceh = R*l and R is the Ricci scalar con- 
structed from the Riemann curvature tensor and *1 is the invariant volume element. The action for non-gravitational 
fields coupling to gravity is denoted by It turns out that any theory which claims to generalize or modify Einstein's 
GR must include classically well-tested GR as a special case and reduce to GR (and hence Newtonian gravity) in the 
weak gravity regime [20] . 

In our work, we shall not address problems arising from certain modified gravitational Lagrangians or their coupling 
to matter or the phcnomenological implications of these theories but scrutinize methods of the derivation of field 
equations that follow from a given modified gravitational action. The dimensionality of spacetime turns out to be a 
crucial parameter for the gravitational models, hence we shall consider some popular three dimensional cases after we 
derive general formulae. 

Relative to an orthonormal coframe, metric variational derivatives are induced by coframe variational derivatives 
and that the vanishing torsion condition can in fact be implemented into the coframe variation either by introducing 
Lagrange multipliers or directly incorporating this condition into the variation by using the variational derivative 
SO" = 0, which relates the variational derivatives of the connection and coframe 1-forms. Using either of the methods, 
one ends up with the same set of metric field equations. However, the latter scheme gets increasingly cumbersome 
with increasing complexity of Lagrangians at hand. Therefore, we shall use outline the Lagrange multiplier method 
appropriate to our investigation. 

In summary, we will work out the Riemannian content of the generalized/modified gravitational Lagrangians which 
contain the scalars built out of the curvature tensor and its varioiis contractions. Specifically, we will consider the 
gravitational Lagrangian densities of type C^q ^ = f{S) * 1 where S is any scalar built out of curvature tensor and 
the function / is an arbitrary differentiable function of its argument. In that, we shall assume that n-form <S * 1 is 



II. CONSTRAINED VARIATIONAL DERIVATIVES 




(1) 
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expressible in terms differential forms such as coframe and connection 1-forms, curvature 2-forms and their contraction 
using Hodge dual operator. As will be apparent, for practical convenience, we adopt the Lagrange multiplier method 
[21-26] using differential forms language (summarized in the appendix for convenience) adapted from [27] throughout 
the paper. Although the Lagrange multiplier method can also be carried out relative to a coordinate coframe, the 
use of differential forms enhances the practical use of the Lagrange multiplier method. Using the method that we 
adopted to derive metric field equations for f{R) theory. Palatini f{R) equations can be derived from the same set of 
variational derivatives, which can be considered as one of the advantages that makes the Lagrange multiplier is more 
efficient method. We shall work in n > 3 dimensions, and a gravitational action Iq has the general form 



Ig = I - 



= [ f{S) * 1 (2) 
Ju 



(n) 

in a local chart U CZ M of n dimensional pseudo-Riemannian manifold where Cq is Lagrangian density n-form which 
already includes the invariant volume form and iS is a scalar which can be expressed in terms of curvature tensor. The 
metric field equations that follow are chart independent provided that the basis coframe and the connection 1-forms 
have the following transformation properties 



H^r = A«^6"', (3) 

-l\a 



0.", ^ (I;", = (A-i)" c^.A'^, + (A-i)" ciA^,, (4) 



for overlapping charts, where the functions A"^ take the values in the Lorentz group. In our presentation, we shall 
suppress the integral sign and use Lagrangian density n-form for convenience. We start with deriving a general 
equation for the variational derivatives of generalized gravitational actions. The basic idea simply is to use product 
rule for the variational derivatives of the scalar fields and subsequently convert the variational derivatives of the scalars 
into variational derivatives of expressions involving forms which allows us to adopt orthonormal coframe. Relative 
to an orthonormal coframe, the total variational derivative of the n— form L^q^ can be calculated as follows. The 
Leibnitz rule for variational derivatives of the product of two scalars gives 

5{f * 1) = {5f) * 1 + /(5 * 1, (5) 

where 5 denotes variational derivative. The variational derivative of the function / can be written as 5f = f'6S where 
/' = ^ . Consequently, the variation can be written as 

SC'^^^ = f'{5S)*l + f5*l. (6) 
Now, using the product rule for the variation 

5{S * 1) {6S) * 1 + 5<5 * 1 (7) 

for the first term on the right hand side of Eqn. (7), Eqn. (6) can be written as 

6£P = f'5{S*l) + {f-f'S)5*l. (8) 

Note here that the second term on tht; right hand side of Eqn. (8) involves only the coframe variational terms whereas 
the variational derivatives of the gravitational variables, such as the connection 1-forms, come from the first term. 
The second term on the right hand side of Eqn. (8) can be simplified by recalling that relative to an orthonormal 
coframes 6°', for a = 0, 1, 2 ■ • ■ (n — 1), the volume n-form can be written as 

* 1 — aaia2a3—an (q\ 

where £aia2a3 - a„ Completely antisymmetric permutation symbol. Using the expression (9), the variational derivative 
of the volume form takes the form 



Note that, unlike exterior derivative, the variational derivative 6 does not pick up a minus sign passing through a 
coframe 1-form and that 5 does not commute with *. Thus, relabeling the dummy indices a\ in (10) and the definition 
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of Hodge dualilty for basis coframe 1-forms, the expression for the coframe variation of the volume element in Eqn. 
(10) can be reduced to 

5*1 = 59°' ^*e„. (11) 

Finally, using Eqn. (11) and Eqn. (8), the expression for the total variation of the action density takes the 
following convenient form 

^4"^ = f5{S * 1) + 56'^ A if- f'S) * 9a. (12) 

This key equation can be regarded as a kind of master formula which has a wide generality, as will be shown below. 
It expresses the variational derivative of the expression f{S) * 1 in terms of the simpler variational derivative n-form 
<S * 1 plus additional terms. For instance, for the simplest case S = R, which will be discussed in the next section, 
in which the Lagrangian density becomes = f{R) * 1 and the total variation will given by Eqn. (20), is the 

simplest example of the general expression (12). It is possible to consider slightly more general Lagrangian density 

where / is a function of a number scalars, say for example, Cq^ = f{S,TZ) * 1. In this case, the total variational 
derivative takes the form 

^4") = ^5{S * 1) + ^5{n * 1) - 59- A (1^5 + - /) * 9a, (13) 

where / is now assumed to be an analytical function of two distinct curvature scalars S and TZ. Not to burden the 
presentation with unnecessary complications, we shall mainly consider the gravitational Lagrangian of type /(<S)*1. In 
order to evaluate the total variational derivative of (12) with respect to independent connection and coframe l-forms 
relative to an orthonormal basis explicitly, the first term on the right hand side of (12) is to be expressed in terms of 
coframe 1-forms, connection 1-forms, curvature 2-form and contractions of curvature 2-forms. The second term on the 
right hand side of (12) directly contributes to the coframe variational derivative and hence to the metric field equations. 
In the Lagrange multiplier method we shall use, the independent field variables arc the coframe and connection 1-forms 
as well as the Lagrange multiplier (n — 2)-forms. However, the connection 1-forms are constrained to be torsion free 
using a Lagrange multiplier (n — 2) form in n-dimensional space. The vanishing of the the torsion cannot be imposed 
straightforwardly to the variation derivatives since it is a dynamical constraint. In such a constrained theory, the 
coframe variations get contributions involving Lagrange multiplier forms. 

The equivalence of the Lagrange multiplier method with the metric variation relative to a coordinate coframe has 
been shown in [28]. However, it should be emphasized that the Palatini method and the metric method, except for 

(n) 

the only the case of the Einstein-Hilbert action Ceh^ lead to different field equations. 

As remarked before, relative to an orthonormal coframe metric equations are induced by coframe variations and the 
equations obtained by independent connection variation can be used to solve the Lagrange multiplier {n — 2)-forms 
in terms of other fields. The Lagrange multiplier n-form terms, which contains vector-valued Lagrange multiplier 
(n — 2)-forms A", of the form 

£LM[A„,w"^,r] = A e« = A (dr + a 9^), (14) 

and Xa is promoted to a new field variable. The variation of the total action is C[g} = Cq\,jj + Clm with respect to 
the independent variables 9°^, uiap and Xa- The metric field equations of the original Lagrangian corresponds to the 
metric field equations of a section of the Ctot determined by 0" = and the field equations obtained from connection 
variation allows one to express Xa in terms of the other fields. The total variation of the total Lagrangian density 
with respect to the independent variables with the torsion constrained to vanish has the general expression 

54":} = 59- A + DXa) + 5ua0 A {n«^ - i(r A A'^ - A A«)} + 5Xa A 6", (15) 

up to an omitted total exterior derivative and the antisymmetric {n — l)-form H"'' = 5Cq^ /6u)a/3 is defined for 
convenience. The metricity condition for the torsion free connection relative to an orthonormal coframe can be 
expressed in terms of connection 1-forms as cJq,^ -|- ujpa = 0. Therefore, the metricity constraint can be implemented 
by anti-symmetrizing the coefiicient of daja/s in expression for the total variation. However, a nonmetricity constraint 
for connection can also be imposed by introducing another Lagrange multiplier (n — l)-form = p^a with the term 
p"^ A (wa/3 + W/3a) = p"^ /^Drjaff added to the original Lagrangian density. The Lagrange multipliers pa0 can be found 
from the symmetric part of the variational derivative with respect to the connection 1-form in this case. We will not 



5 



investigate the modified theories of gravity with torsion and non-metricity (metric-affine gauge theories of gravity) 
[22]. In this work, we simply assume that the connection is metric compatible and therefore it has the anti-symmetry 
property, oja^ + ujpa = and accordingly SuJap + Sujfja = in all the total variations we consider. Consequently, 
one ends up with the property liap + n^a = 0. There is an equivalence between the antisymmetric liap which are 
tensor- valued (n — l)-forms and the Lagrange multiplier A" which are vector- valued (n — l)-forms [22]. The field 
equations SC^^/SiOa^ = can uniquely be solved for the Lagrange multiplier (n — 2) forms as 

A" = 2ipnl''^ + A (16) 

by calculating two successive contractions of the equations SC[g} /Scuais = 0. For n = 3, the wedge product in the 
second term (16) is unnecessary since A" is 1-form in three dimensions. The resulting expression then can be used in 

the DX" term for the metric equation 5C^^/56a = 0. The trace of this term, subject to the constraint 9" = 0, has 
the following general expression 

Oc, A DX" = 2di0{ea A n«'5), (17) 

which applies to all the cases we consider. The equations SC^l/dX^ = O** = DO'^ = imply that the resulting 
covariant exterior derivative D is torsion-free and metric compatible, and therefore D^9" = fi"^ A9^ = 0. 

In the following sections, we shall explicitly work out the Riemannian content of the field equations obtained from 
metric variational principle for various generalized/modified gravitational actions and we shall assume that zero-torsion 
and metricity conditions are imposed for all the models. In particular, we make use of the general expressions (12) and 
(15)-(17). On the other hand we shall indicate how the field equations can be obtained adopting the Palatini approach 
in comparison. If a matter field couples not only to metric components [29] but also to the connection, the resulting 
theory is known as Riemann-Cartan theory where torsion is produced by spin density of the matter [25]. However, 
as has been noted before, torsion may also be generated by gradient of scalar field of a scalar-tensor theory, such as 
BD theory in which BD scalar field couples non-minimally to the curvature. On the basis of equivalence between 
f{R) theory and the BD type scalar-tensor theory one also expects this to be the case for generic Palatini f{R) 
theories. Although we shall not address the problem of formulating Riemann-Cartan type theory for the gravitational 
Lagrangians coupled with spinor matter Lagrangians, our unified approach is even appropriate to accommodate such 
an extension. 



III. SOME PREVIOUS RESULTS: f{R) THEORY 

As an illustration of the general method of calculation for variational derivative of Lagrangian densities that we 
shall present, we first derive the field equations for the well-known f{R) theory. In this section, we shall study both 
metric and Palatini f{R) theories in relation to respective scalar-tensor theories. 



A. Metric f{R) Theory 

We start by recalling that, relative to an orthonormal coframe, the Einstein-Hilbert n-form Lagrangian density can 
be written in the following equivalent forms 

4"l = 5f^"/3A*r'5 = ii?*i, (18) 

where *1 = ffOi---(n-i) jg ^j^g oriented volume element, R = ivi^i^^'^ — ivR^ is the scalar curvature and the dimension 
of the spacetime is assumed to have the values n ^ 3. From purely mathematical point of view, the Lagrangian 
density (18) can be assumed to have the generalization of the form 

4"k = /(^)*l> (19) 

where the function / is assumed to be an arbitrary algebraic scalar function that is differentiable with respect to its 
argument. Similar to the formulae given in (12), the variational derivative of the total Lagrangian density SCq'1;jj 
can be written as 

S^gIh = f'SC'il - 66" A {f'R - f) * (20) 
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where /' = The total variational derivative of the total Lagrangian density £tot i '^"^ j ^a] = ^geh ^lm with 
respect to the independent variables can be written as 

+ 60^ A {f'n„0 A + (/ - Rf) * r + i^A"} + <5A„ A 9", (21) 

up to a disregarded closed form, and from this expression one can read off 

= £)(/' (22) 

Using the general formula (16) and the result (22), the Lagrange multiplier can be found as 

A" = 2 * {df A 0"). (23) 

Consequently, using this expression for the Lagrange multiplier in the equation obtained from the coframe variations, 
one finds the metric field equations induced by the coframe variations as 

- /' * G" + i(/ - Rf) A *r + i:» * {df A r ) = o, (24) 

which was actually derived long ago [14]. Here *G" = —^^^u A *9°'^*'' is Einstein {n — 1)— form. An equation for the 
scalar function f{R) can be found by tracing the metric equation (24). This gives 

(n-l)d*d/'= (^/-/'i?) *1, (25) 

which is identically satisfied for /(i?) = R for the vacuum case. Moreover, for the Einstein-Hilbert case where 
f{R) = R, the Lagrange multiplier vanishes identically since A" = 0. For this case, the Palatini formalism leads to 
the same metric field equations as in the metric formalism (see Eqns. (27) and (28)) as well. We also note the special 
case f{R) = R^ for which the trace equation (25) reduces to the homogeneous wave equation for scalar curvature 
d * dR = in four spacetime dimensions n = 4. However, i?^ * 1 term survive for n = 3 and for n > 4. As can be 
observed from the field equations presented below, the spacetime dimension is usually a crucial parameter that may 
determine physical features of the model at hand. Note that, as expected on consistency grounds, for f{R) = R, 
Eqns. (24) yield the Einstein vacuum field equations. However, unlike Einstein field equations, they contain 4*'' order 
partial derivatives of the metric components relative to the coordinate coframe because of the presence of the term 
D * {df A d'^). This term drops in the constant curvature case il'^'^ = kO^^'^ In this case, scalar curvature is constant, 
R = kn{n — 1) and df = identically. Thus, it is possible to find a Lagrangian density f{R) for which isotropic spaces 
are vacuum solutions. Inserting Q'^'^ = kd'^'^ into (24), one finds f{R) = e'^f"^^)-'^. Cosmological dynamics based on 
the exponential gravitational Lagrangian is studied in [30] and in generic f{R) models R = Q does not obviously imply 
i? = as in GR. In particular, we note here that, for another simple case f{R) = R^ which we shall refer to in the 
succeeding sections, Eqn. (24) yields 

-R^G" - ~ R^ * r + i(TO - 1)(to - 2)di? A r * di? + i(m - !)£> * {dR A 0") = 0, (26) 

which holds for n ^ 3 spacetime dimensions. The solution to the Einstein vacuum field equations R°'^ = trivially 
solves the field equations derived from Lagrangian density f{R) = i?™. Thus, for example, static spherically symmetric 
solutions to the Einstein vacuum field equations are also solutions to (26). But the converse is not necessarily true. 
Not all the vacuum solutions to f{R) theory are also solutions to the vacuum Einstein field equations. 

B. Palatini/Riemann-Cartan Type f{R) Theory 

Riemann-Cartan type gravitational theories accommodate torsion, which is assumed to be generated by spin [31]. 
However, in addition to spin, scalar-tensor theories also give rise to torsion in terms of gradient of a scalar field [32] . As 
we shall show below, for generic f{R) theory, the theory admits an algebraic torsion and consequently the Riemann- 
Cartan connection can be expressed in terms of Riemannian connection and a contorsion 1-form. In this case, the 
metric field equations can be cast into a form completely expressed in terms of Riemannian geometrical quantities. In 
the Palatini method to calculate the variational derivative to f{R) theory, the field equations for connection is used 
to solve the connection in terms of Riemannian connection and additional terms resulting from algebraic torsion in 
terms of the gradient of the scalar df . The Palatini approach to f{R) theory, which we shall study in this section, is 
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therefore can be considered as modified gravitational theory with algebraic torsion. We shall start with formally the 
same Lagrangian density (18). In the case where non-zero torsion is allowed [33, 35], the field equations, which can 
be found from the variational derivative (21), are 

df A *e°'^ + f'e" A *e"^i, = o, (27) 

f'*G" + {f- f'R) *e°'=0, (28) 

respectively. Note here that for f{R) = R, the Palatini variation yields 11"'^ = identically if one imposes ©" = 
after the variation in the case of C^^h ■ This follows from the fact that in the metric variation for the Einstein-Hilbert 

Lagrangian, one finds from (23) that A" = [36]. As a result Palatini and metric variational derivatives lead to 
the identical metric equations in this exceptional case of the Einstein-Hilbert Lagrangian. The field equation which 
follows from the independent connection variation, Eqn. (27) is an algebraic equation for the torsion 2-form ©" which 
can be solved for it as 

e'' = j^O'^^dlnf, (29) 

using the identity 6" A iaX = pX which holds for arbitrary p-form A. This result then can be used to solve the 
corresponding connection 1-forms which can be obtained from the Cartan's first structure equations. 

Returning to the Palatini equations, the variational derivative with respect to connection, namely H"^ = lead 
to an algebraic torsion in terms of the gradient df . As a result, in Palatini formalism, even if the connection is 
metric compatible (and therefore commutes with contractions), it has torsion. However, algebraic torsion allows 
the elimination of the independent connection 1-form and the field equations can be cast into a metric theory with 
additional gravitational currents under suitable conditions stated below. For this purpose, it is convenient to introduce 
contorsion 1-form = —Kp^ defined by 0" = K"'^ A 9^^ . The contorsion 1-forms can be expressed in term of the 
exterior derivatives of the function /' as 

K"^ = (;i32) {O'^idlnf')^ - 0^(dln/')"} , (30) 

where the indices of the 1-form d In /' are relative to the orthonormal coframe we employ. In terms of contorsion 
1-form, the metric compatible connection 1-form w"^ can be written as sum of Riemannian connection co"^ and a 
contorsion term as 



u''^=w"0+K"^, (31) 

where, in order to distinguish metric connection from the Riemannian connection, the dotted symbols for the Rie- 
mannian quantities are used. The result (31) will be elaborated below. 

The Cartan's second structure equations for the corresponding curvature 2-form then lead to similar decomposition 
of the curvature 2-forms as well. One finds 

^ Qa0 ^ £)^a/3 ^ j^a^ ^ ^m/?^ (32) 

where is the Riemannian curvature 2-form and D is the covariant exterior derivative corresponding to Riemannian 
covariant derivative. Note that / and hence /' is a function of scalar curvature R. Thus, the elimination of independent 

connection requires R to be expressed in terms of the dotted quantities using the trace of the metric equations (28). 
If we assume that this requirement can be satisfied, then the Einstein (n — l)-form can be written as 

f*G"=f*G"-D* {df A r) - 1^^-^ * T^if], (33) 

where the energy momentum forms for the scalar field /', namely 

* Taif] = -\ {{iadf) * df + df A i„ * df} , (34) 

is defined for convenience. Note that this result explicitly shows that non-zero torsion produces the kinetic term for 
the function /' at the level of action and this corresponds to the effective energy momentum term *T^[f] in the 
resulting metric equation. Finally, using this result and dropping the dots over Riemannian quantities, the metric 
field equations can be written as 

/' * G« = i(/ - Rf) *6" + D* {df A r) + J^^^ * T«[/'], (35) 
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which arc expressed in terms of Riemannian connection 1-form and the corresponding Riemannian curvature tensor. 
Taking the trace of the metric field equation one finds 

{n-l)d*df' = ^{n-l)ydf' A*df' + yR*l (36) 

Note that because of the term D * {df A 9"), the field equations (35) are fourth order partial differential equations 
relative to a coordinate basis. The above reduction of f{R) theory with torsion to a metric theory has been carried 
out recently using the equivalence of such a theory with Scalar- Tensor theories relative to a coordinate basis and it is 
also argued that the theories with torsion and/or non-metricity can also be cast into metric theory for f{R) theories 
[33]. As we shall show explicitly in the following sections, this is not true for other modified gravitational Lagrangians. 

We now briefly discuss the equivalence of generic f{R) theory with torsion to the Scalar-tensor theory. Note that 
if the field variable f'{R) is promoted to a scalar field, then equations of motion (24) and (25) for f{R) theory have 
formal resemblance to the metric equations in BD theory except for energy momentum (n — l)-form corresponding to 
the kinetic term for the BD scalar function. In other words, the equations of motion for generic f{R) theory can be 
obtained from a Brans-Dicke type action with kinetic term for BD scalar that minimally couples to the metric tensor 
is absent. Consider the following BD-type Lagrangian density in n dimensions 

C = Cbd + V{(1>) * 1 = A *e^"' + ^d(j) A *d(f> + VU) * 1, (37) 

2 2(p 

where the scalar field (j) couples to the metric non-minimally and there is a self interaction potential term for the 
scalar field. The metric field equations, that follow from (37), is 

* G" = * T"[</)] +D*{d(pA r ) + y((/>) * r, (38) 

where the canonical energy momentum (n— l)-form field *T°'[<j)], which results from commuting the coframe variation 
with Hodge dual, is 

* Ta[4>] = -|{(ia#) * d<i) + d(l) A ioc * #} (39) 
whereas the field equation for scalar field can be written as 

{^ + ^)d*d(l) = V' *1. (40) 

For the BD theory with vanishing BD-parameter, namely for w = 0, the metric and the scalar field equations that 
follow from the variation of the action for this Lagrangian with zero torsion condition imposed can be written as 

-<I)*G" +D*{d(pA0'') + V{<p)*9" =0, (41) 
R = V'icP), (42) 

respectively [34]. Here, ' denotes derivative with respect to the argument of the scalar potential, V = The trace 
of the metric equation (41) and (42) can be combined to have 

(n - l)d * d(f> = {nV{(p) - \(n - 2)(1>V\4>)) * 1. (43) 

Comparison with the metric field equations for the /(i?) theory, assuming that / is a definite functional form and 
that /" 7^ 0, the equation f{R) = (j> can be inverted in order to be express R in terms of cj). Then, R{(j)) can be used 
in the equation 

nV{<i>) - i(n - 2)cl>V'{4>) = f / - f'R, (44) 

to determine the scalar potential V{(j)). In particular, for n = 4 one has 1^((^) = \ f{R{(j))). As a result, with the above 
identification and redefinition of the field variables, the field equations for f{R) theory transformed into a BD-type 
theory. In the BD-type scalar theory with field Eqns. (41) and (42), the energy momentum (n — l)-forms for any 
matter field *T'^[$] that couples only to the metric field but not to the derivatives of the metric field, are covariantly 
constant. Explicitly, using the contracted second Bianchi identity, which is a purely geometrical identity, namely 
D * G" = 0, the covariant derivative of (41) gives 

- # A *G" + A*{d(pA0^) + V'{4>)d(p A *e" = 0, (45) 
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by virtue of the field equation (42). In other words, D * T^[$] = fohows from the fact that the covariant exterior 
derivative of the vacuum metric field equations (41) vanish. Conversely, the metric field equations (41) yield the 
field equation (43) by taking the covariant exterior derivative of (41) and using (42). The fact that the field equa- 
tions are covariantly constant can be shown to follow from more general property of Lagrangian density that it has 
diffeomorphism invariance, independent of the resulting field equations. 

It is important to emphasize again that, the Palatini and the metric variational principles in general, lead to different 
sets of metric field equations. The scalar tensor equivalence for the corresponding f{R) theories above provides an 
example for the inequivalence of the two theories since the Palatini variation f{R) theory lead BD-type scalar tensor 
theory with BD parameter u) = — ("I2) with potential term for the scalar field whereas f{R) theory based on metric 
variational principle lead to scalar-tensor theory with cj = [33, 40]. 

The Palatini connection is conformally related to Levi-Civita as in (31) with the contorsion form expressed in 
terms of a scalar gradient (30). Because the conformal factor involves scalar curvature, the parallel transport of the 
connection in the Palatini method leads to a distinct autoparallel curves which can be expressed in terms of Levi- 
Civita connection. A timelike vector field X, which is normalized as g{X,X) = —1, is tangent to autoparallel curves 
of the connection (31) if it satisfies 

WxX = 0. (46) 

Here, it is convenient to introduce the operator b which maps the vectors to their metric duals, i.e. to the corresponding 
1-forms. It can be defined implicitly in terms of the metric tensor as X^{—) — g{X, — ). Since is metric compatible 
and hence commutes with the operator b, we have VxX^ = {VxXy. Using the expression (31) for the connection 
1-form, Eqn. (46) can be written in terms of the Riemannian covariant derivative Wx as 

VxX^ + ^ix{dRAX^)=0. (47) 

As a result, similar to the scalar tensor theories [37], the non- Riemannian connection given in (31) induces an 
acceleration term depending on the gradient of the scalar curvature into the autoparallel curves. We note here that 
the geodesies are determined by the Riemannian structure instead of the independent Palatini connection as a result 
of the fact that the matter energy momentum (n — l)-forms are covariantly constant: D * = [38]. On the 
other hand, the correct cosmological behavior of /(i?) theories in the Palatini formalism requires consistent averaging 
procedure of energy momentum forms describing the microscopic(atomic) structure of matter fields and it also predicts 
alterations to atomic structure and particle physics laws, see e.g. [39]. 



C. Equivalence to Scalar-tensor Theories at the Level of Action 

The equivalence of f{R) theory to particular Scalar-tensor theories discussed above by comparing the corresponding 
field equations can also be deduced, perhaps in an indirect way, at the level of action. (Such an equivalence holds for 
f{R) theory with/ without torsion.) The equivalence can be achieved by introducing an auxiliary scalar field x and a 
scalar Lagrange multiplier A imposing the constraint x — = to the gravitational action. 

Wc shall mainly consider modified gravitational Lagrangians of the form f{S) * 1, where / is an algebraic function 
of scalar <S and we assume that 5*1 can be expressible in terms of basis coframe forms, connection 1-forms, curvature 
2-forms and their various contractions. In the place of f{S) * 1, consider an auxiliary Lagrangian density Caux = 
J0.aux[x, \ of the form 

jCaux = [f{x) + Kx-S)]*l. (48) 

The total variational derivative of the auxiliary Lagrangian (48) with respect to the independent variables indicated 
above yields 

SCaux = [Sxif + A) + ^A(x - 5)] * 1 - X6{S * 1) + SO'' A (/ + XS) * 6^ (49) 

where S denotes the variational derivative as before and the prime denotes derivative with respect to the auxiliary 
field X- The equation for the scalar field x immediately yield the Lagrange multiplier A = — /'. Replacing this result 
and the corresponding constraint back into the action (48) one finds the dynamically equivalent Lagrangian density 
as 



Caux-egv. = [/(x) + (X " 'S)/'(x)] * 1, 



(50) 
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up to an overall sign. The resulting Lagrangian density Caux-eqv. now has no dependence on the Lagrange multiplier 
A and in the equivalent Lagrangian density /' couples to the scalar S. However, the reduced form of the auxiliary 
Lagrangian density allows redefinition of the field variables. For f{R) theory, namely for S = R, Eqns. (49) and (50) 
reduce to 

SC^ux = [Sxif + A) + SX{x -R)]*l- \S{R * 1) - 5w^^ A D{X * 6^ 

+ A {-Xn^. A + (/ + XR) * 0^} , (51) 

and 

Caux-egv. = f'^a0 A ^^"^ + (/ - f'R) * 1, (52) 

respectively. Although the auxiliary action (52) consistently simplifies to the original action, this form of the action 
further suggests the redefinition /'(i?) = </>. This allows one to bring the action (52) into equivalent Scalar-tensor 
type action given in (37) with w = and with an appropriate potential term for the scalar field assuming that (j) can 
be written in terms of R. This general mechanism of generating non-minimally coupled scalar field in generalizing 
gravitational actions will be used to propose a modified Chern-Simons theory in four dimensions in a later section. 

IV. VARIATIONAL DERIVATIVES OF MODIFIED GRAVITATIONAL ACTIONS 

In addition to the generic f{R) theories, modified gravitational theories based on the Lagrangian densities involving 
quadratic curvature invariants and their algebraic generalizations also lead to fourth order models. These models have 
a long history [41] , they have been studied in different contemporary contexts ranging from cosmology and perturbative 
approach to quantum theory gravity [42]. This section is devoted to the calculation of the metric field equations for 
fourth and sixth order theories. 



A. Generalized Gauss-Bonnet Actions 

The outcome of employing extra dimensions for unification as well as quantum corrections of gravity requires some 
well behaving generalizations of the EH Lagrangian, such as scalars constructed from of higher order curvature terms. 

Actually, using an arbitrary scalar in action constriicted from Riemann tensor leads to pathologies [43, 44] such as 
to have more than required degrees of freedom in equations of motion or existence of ghosts. However, there is a 
special combination of such terms in action, as shown by Lovelock [45], which are known as dimensionally continued 
Euler forms or Lovelock scalars. It turns out that, the Lagrangians constructed from these terms have some desirable 
properties such as, though nonlinearly, the field equations involve at most second derivatives of metric components 
[45, 46], theory is free of ghosts [43, 48], i.e. kinetic terms in the action with the wrong sign, and in four dimensions 
the only surviving term apart from cosmological constant is the EH Lagrangian [45, 49]. This generalization of 
considering Euler forms were studied in different contexts [50-55] such as Kaluza-Klein reduction, low energy limit of 
string theory, brane-worlds, and highcir dimcinsional exact sohitions. 

Thus, in higher dimensions, one can consider, in addition to the usual Einstein-Hilbert term, the combination of 
Euler forms 

[«/2] 

7g=^4''"\ (53) 

fe=0 

where 

4'''"' = ^c.,h A • • • A fia./3. A (54) 

For k — 2n, is an exact form and is proportional to Euler- Poincare characteristic of even dimensional manifold. 

For simplicity, we consider the gravitational Lagrangian density n-form 

4") = A fl^, A 

= 2fi„/3 A - ARo, A + i?2 * 1, (55) 

corresponding to fc = 2 for n ^ 4. This term is called as Gauss-Bonnet or Lanczos Lagrangian, it is a pure divergence in 
d = 4, and first introduced by Lanczos [46] . Although here we only consider Gauss-Bonnet term, it is straightforward 
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to extend the formalism to the higher order Euler form Lagrangian densities. Now. it is convenient to define the scalar 

G = RafffiiyR""^^^ — '^RfxvR'^'' + R? which separates the volume form so that C^g ^ = Q *1. In four dimensions, £g ^ is 
a total differential and does not contribute to the field equations. This follows from the identity 

^e^^^flf"" = Q"!^ - i (61" A i?'^ - 61^ A R") + {RB'^f^ , (56) 

which relates the left dual of the curvature 2-form to the curvature 2-form itself and its contractions [47]. From (56), 
it follows that 

^lc.p^n^^,^*e""^"' = -\n''0^Q.|^^, (57) 

as a result of first Bianchi identity and the right hand side is an exact form (see Eqn. (107)). Therefore, the GB term 
of the form 17^/3 A A *6"^'^'' does not contribute to the field equations. However, it is possible to make C'g^ to 
contribute to field equations, for example, by coupling it with a scalar field. On the other hand, the presentation of 
f{R) type Lagrangian densities indicates that f{Q) * 1 type Lagrangian density will also have similar features that 
/' will act like a non-minimally coupled scalar field. After these preliminary remarks, we will derive the equations of 
motion for the modified Lagrangian density of the form 

CPg = f{g)*l, (58) 

for n > 4. The independent variations with respect to the coframe and connection 1-forms axe constrained to be 
torsion free by Lagrange multiplier term (14). Using the same steps in the derivation of the variation (12) from the 
modified Lagrangian (19), the total variation of the Lagrangian density, 

4:t^=41+'CiM, (59) 

is related to the variational derivative of the Lagrangian density £g"' as 

SC[:l = /'(5£^") -S0'^A{fg-f)*0^ + 6£LM. (60) 

Note that for n = 4, there is no contribution from the term containing SC^g'^ on the RHS of (60). Thus, the cases 
n > 4 and n = 4 lead to different field equations. The variations on the right hand side of this expression can be 
calculated to find 

5L^t^l = Scoafs A {n"'^ - i(r A - 6"^ A A")} + <5A„ A 6" 

+ Se^ A A A P * 9"""^ + if'G -f)*e^ + DX^}. (61) 

Here ' denotes derivatives with respect to the scalar Q and the expression for Hq,/3 obtained from the variational 
derivative with respect to the connection 1-form yields 

n"'^ = 2D (/'fl^, A ^e^""^) , (62) 

which is to be evaluated subject to the vanishing torsion condition. In contrast to the other field equations that we 
shall consider below, the field equations for the GB terms are second order in metric components. This follows from 
the fact that the second Bianchi identity satisfied by the curvature 2-forms and the Lagrange multiplier terms can 
be expressed in terms of curvature 2-forms. The order of the metric equations persist even in the corresponding first 
order formalism where an independent dynamical torsion is present. 

As a result, the expressions for the Lagrange multipliers can be inserted into the metric equations 5C^^l/59x = 
to find 

fn^, A A *9'^''^p^ + (/ - f'G) * - a d * {df a 0^^^^) = o, (63) 

for n > 4 dimensions. The field equations (63) that follow from the f{G) Lagrangian axe similar to the field equations 
for the Lagrangians with a scalar field which couple non-minimally to the GB-terms. The term Oj^j, AD* {df A 9'^"^) 
contains 4*'* order partial derivatives of the; metric components relative to a coordinate coframe. Now specializing to 
n = 4 dimensions, in particular, to the gravitational Lagrangian c['^^ = £^]j + jCqq^j the metric field equations can 
be written as 

- *G" + (/ - f'G) * 9°' - A^^^ A D * {df A 6""^") = 0. (64) 
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Apparently, the order of the metric field equations for f(Q) augmented to four as a result of algebraic generalization 
of GB term. We refer to the review article [17] for f{Q models and their applications. 

Before we begin to study the the individual terms appearing in (55), we note that the dimension of spacetime is 
crucial and in particular, in four spacetime dimensions, using (55) and (56) the most general Lagrangian which are 
quadratic in the scalars the curvature can be written in the form 

=ai?" A*i?a + 6i?2*l, (65) 

with a, b being constants. The Lagrangian density (65) is the most general gravitational Lagrangian density involving 
quadratic curvature terms in three spacetime dimensions as well and we will study of this special case in a later section 
below. However, the most general Lagrangian involving quadratic curvature terms is to include term in 

n > 4 dimensions. 



B. Ricci-Squared Action 

In this subsection we will calculate the variational derivatives of the Ricci-squared Lagrangian density whereby the 
corresponding field equations for the gravitational Lagrangian densities in (65) with arbitrary a, b and its algebraic 
generalizations of the form f{R,V) can be found using the results of the previous sections. We first consider a 
Lagrangian density containing Ricci tensor-squared term, which can be written in terms of Ricci 1-forms as follows 

r^,"^ = i?„ A *R" = RapR"'^ * 1. (66) 

(n) 

The metric field equations corresponding to C-p will also be used later in the study of three dimensional massive 
gravity models. We first derive field equations for the Lagrangian density in a form appropriate to our presentation 

and then derive the field equations for its generalization of the form C^Q-j^ = f{V) * 1 where V = Ra/sR'^^ ■ We first note 
that the symmetry of the Ricci tensor indices R"" A9a = follows from the contraction of the first Bianchi identity, 
namely DQ" = fl"p A 9^ = together with the anti-symmetry property il^^ + il^^ = of the curvature 2-forms 
relative to an orthonormal coframe. Thus, the symmetry of the components of the Ricci tensor can be incorporated 
conveniently into the total variation without introducing an additional constraint. We constrain the variation to be 
torsion- free so that we have the total Lagrangian density = + Clm- This Lagrangian has the total variation 

lI"^} = Sua^p A {D * (61" A i?'^ - 61^ A i?") - 1(6*" A A^^ - 6"^ A A")} + A 6" 

+ (56"" A {n^f} A V * (6*" AR^ -e'^ A i?") - v(i?„ A *i?") + DA^} , (67) 

where the identity i?^^ * 1 = A *6'"^ has been used to obtain (67). The resulting vacuum metric field equations, 
subject to the condition that 0" = 0, can be conveniently written in terms of 1-form 7-"" = V^pO^ as *7-*" = where 
*V°' is given by 

^V<' = Ail"* {9" AR^ -9^ A R") - i^iR^ A *R") + 2Di„D * {9" A i?'' - 6"^ A i?"). (68) 

To our knowledge, the metric field equation for the Ricci-sqiiared model formulated in terms of differential forms (68) 
is new. The covariant exterior derivatives in the Lagrange multiplier terms explicitly show that the Ricci-squared 
model is a fourth order model. 1-form defined by (68) is symmetric with respect to its indices Pafi = V^a- We 
also note here that the trace of the metric field equations (68) can be found as 

V\*l = {n- 4)i?" A + f d * dR. (69) 

The trace formula (69) again reflects the fact that the spacetime dimension is important for the Ricci-squared models 
as we shall see in three dimensions. This form of the field equations allows one to find a nontrivial solution to the 
A- vacuum field equations almost by inspection of the metric field equations. For example, i?" = k9'^ with the constant 
k solves the A— vacuum field equations *V°' = —2(n — 1)A;^ * 9". 

Now, having the total variational derivative of £p ^ at our disposal, the field equations for = ./(^) * 1 can be 
found. As before, we also assume that / is an algebraic and differentiable function of its argument. The corresponding 
metric field equations can be found by straightforward application of (12). The total variation of C^^^ = Lq^+SCj^j^ 
then becomes 

6Ci:} = f'6C^^^ - 69^ A {f'P -f)*9^ + 5C^^. (70) 
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Carrying out the total variation using (12) one can find 

56^t = Scoap A {D * f'{R" Ae^ -R^ A 6'^) - \{e" AX^ -6^ A A«)} + A 6" 

+ A {2f'n^0 A V * {R" A e^) - f'i^,{Ra A + {f'v -f)*e^ + dx^} . (n) 

For the purpose of comparison, it is convenient to write the resulting metric field equation using 1-form Va defined 
in (68). This gives 

f'*V°' + DifX") + {f'v - /) * r = 0. (72) 

In contrast to the metric formalism, if the Palatini formalism is used as in [18], then from (67), the equation for the 
connection can be found as 

U""^ =D* f'{9" AR^ -e^'A R^) = 0, (73) 

where D in this case is non-Riemannian covariant exterior derivative. Eqn. (73) has to be solved for connection 
1-forms and then subsequently replaced into the coframe variation to obtain the metric field equations. Obviously, as 
in the case of f{R) theory, it follows from (73) that the Palatini and metric formalism lead to totally different set of 
equations for the metric field. However, adopting Palatini formalism with an independent and torsion- free connection, 
the recent work on the models based on the Lagrangian n-form of the generic form /("P) * 1 lead to the result that 
such models accommodate metrics of Einstein spaces [56] . 

Returning to our treatment of the models with quadratic ciirvatiirc invariants using the metric formalism, we 
now study another popular fourth order model with the Lagrangian density containing Riemann tensor-squared 
(Krectschmann scalar) terms [57], namely the n-form 

= A *n^^ = ^Rc^fi^.R^^^" * 1, (74) 
and its algebraic generalization is of the form 

4")=/(/C)*l, (75) 

with / being an algebraic function of the Krectschmann scalar JC = Rafifn^R"^^" ■ The field equations for this 
Lagrangian, with torsion tensor constrained to vanish, can be obtained from the variation of the total Lagrangian 
density 

4t = 4'^+£LM. (76) 
For the Lagrangian density Cq^, equation (12) becomes 

SC'^l = f'56£^ - 56^ A \{fK -f)*e^ + 5C^j^. (77) 

The total variation with respect to independent coframe and connection 1-forms and Lagrange multiplier (n — 2)-forms 
then can be written as 

56^^ = SuJocp a{D* {ffl"^) + i(A" AG^ -X^ A r )} + 5Xa A 6" 

+ 56^^ A {-!/' (vO„^ A - A V * fi"'') + Wk: -f)*e^ + dx^} . (78) 

The field equations 6C[gl /Sujais = can be solved for the Lagrange multipliers A" using (16) with 11"'' = D{f' * Q,"^). 
Therefore, the metric field equations subject to the condition that 9" = 0, can be written as 

ADi^Dif * n«^) -e^A Dij0D{f' * n"^) + f * T^[n] + \{f'K - /) * = o, (79) 

where the canonical energy momentum tensor for the curvature 2-forms, 

*T^[n] = ^^ (^^0„/3 A *f2"'' - A v * ^'"^) , (80) 

are introduced for convenience. Note that this metric equation holds in n ^ 3 dimensions. For / * 1 = fi^^ A 
and /' = 1, the field equations (79) reduce to 

DiacD * ri"^ - 16*" A Dia^ifiD * f]"'^ + *T^[n\ = 0, (81) 
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which we shall make use of in the following sections. The trace of Eqns. (81) leads to the following equation 

(n - l)d *dR+{n- 4)0"^ A *Qc.0 = 0, (82) 

using the Bianchi identity D * G" = 0. As for Ricci-squared and f{R) = B? models, Eqn. (82) depends of the 
dimensions of spacctimc where n = A appears to be special. 

The method adopted can be used to calculate the field equations easily for more complicated-looking actions 
compared to coordinate methods. As an illustration of this point consider generalized gravitational Lagrangian 
density in curvature scalar, for example, of the form f{Q) = Q"^, 

g'^=K? + 16V^ + R'^-8R'^V- 8K.V + 2]CR^ , (83) 

contains cross terms of the form R?1C, R?1C and KV. The variational derivatives of these types of individual cross 

terms as well as the models based on the generalized Lagrangian of type, for example /(-R, Q), are of some interest. In 
order to illustrate how the method we adopted can easily be used to calculate the variational derivatives of Lagrangian 
densities in curvature polynomials, consider the generalized Lagrangian density of the form 

6^}^ = R^R''^^"'Ra0^.. * 1. (84) 
As before we constrain the variation by adding Clm to the corresponding action and perform a total variation 

<5£^4 = 2i2/C5£^if + R'^SiQa^ A *fl°'^) + SjClm- (85) 

This explicitly gives 

^£^'4 = '^'^"/3 A {D{{ICR - i?2) * r^^} - i(r A A'^ - 61'' A A")} + SK A 6" 

+ Se^^ A {-ICR *Gf,-R^ {if,na0 A - A V * fi"'') + DX^} . (86) 

Solving SC^^2'fi/Sujaf3 = for the Lagrange multipliers A^ in terms of other field variables using (16), and using the 
resulting expression for A^ in the equation SC^^}^^/ 59^ = 0, the metric field equation becomes 

- /Ci? * G" + 2i?2 * T°'[a\ +D* {d{RIC - R^) A 9°'] = 0. (87) 

As for all the generalized gravitational Lagrangians, written in the above form, field equations has formal resemblance 
to the field equations (38) and (40) of BD type Scalar-tensor theory where scalar field is replaced by a function of 
curvature scalars. 



C. Weyl Tensor-squared Action 

Weyl tensor is the trace-free part of the Riemann curvature tensor and therefore the corresponding Weyl 2-form 
Cap satisfy Cap A *9°'^ = 0. Recall that relative to an orthonormal coframe in n > 4 dimensions curvature 2-form 
has the expression [27] 

= - ^ [9- AR^-ef'A R-) + ^^_,l^^_,^ Re^^. (88) 

It is conformally invariant, has the same symmetry properties with the Riemann tensor and exists in spacetime 
dimensions n > 4. There is a unique action based on the contraction of the Weyl 2-form with itself 

C'^^ = Cap A 

= nap A - ^R" A *Ra + * 1. (89) 

In four dimensional case, Lagrangian density (89) is locally scale invariant and the metric field equations derived from 
C^^^ arc know n as Bach- Weyl equations which have 4*^ order partial derivatives of the metric tensor relative to a 

coordinate basis [49]. 

As we have mentioned above, the GB Lagrangian density c'^g^ (55) is a total differential. Thus, without altering 
the equations of motion for , by adding the term Cq^^ to the Lagrangian density c'^) , the term flap A ^ft"^ in the 
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Lagrangian density can be eliminated. This gives somewhat simplified four dimensional gravitational Lagrangian 

density 

/;W = i?« A *Rc - iB? * 1, (90) 

up to a total differential and an overall multiplicative constant. This reduction of the Lagrangian density makes 
it possible to derive the field equations of the model from those of of f{R) = and the Ricci-squared model 
studied above. This requires the use of the properties and the identities satisfied by the curvature 2-form and therefore 

the cases of the three and higher dimensions has to be handled separately. Returning to the general (n ^ 4) case and 
following the preceding examples of generalized Lagrangian densities, it is convenient to define 0-form C = Cn^afiC^""'^ 
and similarly the generalization of Cy^ is assumed to have the form 

4t = /(C)*l, (91) 

where / is an algebraic function of the scalar C. Note that the local scale invariance of the original Lagrangian has 
been lost in such a generalization. We will study the metric field equations with torsion constrained to vanish and 
consider the total Lagrangian density £j"j = J^Q-^y + Clm- The variational derivative, in the light of the variational 
relation given (12), is then given by 

5Ctl = /'(^r +56- A{f-f'C)*e^+ SjClm , (92) 

where ' in this case denotes derivative with respect to the scalar C. Carrying out the variations indicated in this 
equations one finds 

sci^^l = A I :^Rp A + /' * T^[w] + (/ - f'c) * r + dx^ | 

+ SiOocff A {2D{f' * C^) - i(r AX^ -01^ A A")} + (5A„ A 9", (93) 

where the definition of canonical energy momentum tensor for the Weyl 2-form 

* T^[W] = -{iaC^.) A *C^" + C^^ A ic * C"'^ (94) 

has been used for convenience and the irrelevant total exterior derivatives are omitted. The field equations 
dc[^} /()LUai3 = can be solved for the Lagrange multipliers A" using (16) with 11"^ = D{f' * C"^). Therefore, 
the metric field equations subject to the condition that the torsion vanishes, becomes 

f'{j^R0 A *C^" + *T"[W] + ADipD * C^" + 9" A Di^i^D * 

+ (/ - f'C) *e°'+ Difsidf A *C^") = 0. (95) 

As a result, the effect of algebraic generalization of the Weyl tensor-squared action is (i) the overall /' factor for the 
first term and (ii) the last two improvement terms in (95). As expected for consistency of the field equations, for 
/(C) = C the field equations (95) reduce to Bach- Weyl equation [49, 58]. As before, the Palatini variational derivative 
also lead to field equations that accommodate torsional connection for Wcyl-squarcd action [58]. Recently, in n > 4 
dimensions, the corrections to Tangherlini solution [59] arising from higher order curvature terms, namely the terms 
of the form C * 1 are studied in [60] and in this context, it is possible to consider corrections originating from the 
terms of the form /(C) * 1 in four dimensions since in four dimensions the Schwarzschild solution is not modified by 
lowest quadratic curvature corrections of the form i?^ A or i?^ * 1 [60] . 



D. Actions Involving Derivatives of Curvature Tensor 

In this section we briefly discuss the gravitational Lagrangians involving derivatives of curvature invariants of type, 
for example RAR [61], where A = dd^ + d^d is the Laplace-Beltrami operator. The codifferential d^ is the metric 
dual of the exterior derivative d and we shall use the definition d* = {-1)p+^ * acting on p-forms. The models 
based on this type of Lagrangians as well as all of their algebraic generalizations considered above lead to the metric 
field equations which are 6*'' order partial derivatives of metric components relative to a natural coframe. The 
field equations for such theories involving derivatives of scalar curvature terms have found applications in cosmology 
[62, 63]. 
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The calculations below can be extended to the generalized Lagrangians of type dRa0 A*dR°'^ and dRa^^i^ A*dR°'^^'^ 
and also to their polynomial generalizations in a straightforward manner. 

We shall consider the simplest possible gravitational Lagrangian density involving derivatives of the curvature 
tensor, namely RAR [61] and its algebraic generalizations. In order to apply the variational procedure used above for 
generalized gravitational actions presented, it is first convenient to rewrite the action density RAR* 1 in the following 
form 

c'^^^ = yRA*dR=-^{dRf*l, (96) 

equivalent to the expression RAR* 1 up to an omitted total derivative and an overall sign. We define the scalar [dR)"^ 
which takes the form (dR)^ = g°'^ {daR){dpR) relative to a coordinate frame. This definition will be convenient 
in the subsequent generalization of the Lagrangian (96) below. The derivation of higher order field equations for 
gravitational Lagrangians of the form RA'^R * 1 is presented in [63]. 

As before, we constrain the total variation so that the torsion vanishes by adding Clm to the original Lagrangian 
density (96) as 

4ot = A *dR + Clm- (97) 
We found the total variational derivative of the total Lagrangian density to be 

64"! = Sloc^p a {2D((Ai?) * e"^) - i(r AX^ -e^ A A")} + SX^ A 6" 

+ SOa A {{AR)n^^ A ^e^""" - RAR * 6»" + *T"[i?] + DX"} , (98) 

up to a closed form and the canonical energy momentum (n — l)-forms for the scalar curvature R energy momentum 
forms of a scalar field, 

* Ta[R] = -^{{iadR) *dR + dRAia* dR}, (99) 

has been used for convenience. The canonical energy momentum form of the scalar curvature results from the 

commutation of variational derivative of the basis coframe 1-form with the Hodge dual exactly as in the case of scalar 
field (cf. Eqn (39)). The resulting metric field equations [61, 62], subject to the condition that the torsion vanishes, 
are given by 

(AR) *R" + l* T"[R] -D* {{AdR) A r } = 0. (100) 

These equations are 6*'' order in the components of metric tensor since turns out to have terms which are 5*'' order 
in the partial derivatives of the metric components. The trace of the metric field equations is given by 

*{A + R)AR- ^^^dRA*dR^O. (101) 

Both (100) and (101) have also formal resemblance to field equations of a BD type scalar-tensor theory (38) and (40) 
respectively. If matter fields couple minimally to gravitational action (96), then the "scalar field" AR act as a variable 
coupling constant for the metric tensor. In this case, the trace of matter energy momentum fields acts as a source 
term for the wave equation for the scalar AR in (101). 

As in the case for the f{R) theory studied above, for the sixth order theory based on (96), the Palatini variational 
derivative lead to a metric field equation different then (100) as well as it induces a torsion for the connection. In this 
case, the connection 1-form can be found by solving the relation 

D{{AR) * 0"l^) = 0, (102) 

which lead to nonvanishing, algebraic torsion similar to the case for Riemann-Cartan type f{R) theory. 

Finally, we will consider the algebraic generalization of the sixth order theory based on (96). The model then has 
the Lagrangian density of the form 

^gL = fiidRf) * 1, (103) 

where, as in all the cases before, / is assumed to be an algebraic function of the scalar (dR)^. The total variation of 
the total Lagrangian density 



(104) 
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with the 9" = condition imposed by Lagrange multipher term can be found using (12). This gives 

sci:l = f'dc'^^l - se^ A {{dRff - /} * r + sc^^^. (io5) 

Using the total variation (98), the field equations that follow from the expression (105) is 

/'(Ai?) * i?" - i/' *T"[R]-D* {f{AdR) A 6*"} + i {{dRff - /} *e" =0, (106) 

which is similar to the field equations based on the model (96) except for the the last terms in curly brackets and 
the new factor /' as the coefficient of the Ricci 1-form. ' denotes derivative with respect to {dRf. The generalized 
metric equations (106) arc 6*'' order in metric components which is of the same order as the model based on (96). It is 
interesting to note that, in contrast to the fourth order models, the model derived from the generalized gravitational 
Lagrangian density Cq^ = C^eh + * 1 + bdR A *dR is studied in [62] and is shown to be conformally equivalent 
to the conventional General Relativity coupled to two interacting scalar fields. 

V. GRAVITATIONAL CHERN-SIMONS TERM IN FOUR DIMENSIONS 

The generalized gravitational actions studied in previous sections indicate that any algebraic function of curvature 

invariants at the level of action behaves like non-minimally coupled scalar field in the corresponding field equations. 
In this section, in the light of this important observation, we find the field equations for gravitational action based on 
the generalization of gravitational Chern-Simons terms [64-66]. In 4n dimensions the Chern-Simons form in terms of 
1-form-valued matrix notation can be written as 

Tr[f22"] = A n'^^ A • • • A Q''^ A Q-"^ = dK, (107) 

where Tr denotes trace and (4n — l)-form K can be written in terms of connection 1-form as 

K = 2n r dtt^^''-^^Tr[oj{doj + iw2)(2"-i)]^ (108) 
Jo 

as a result of the fact that Tr[a;^"] = [67]. The indices and the wedge products are suppressed for simplicity of 

the notation. Wc specialize to four dimensions, namely n = 1, and in particular to the Einstein-Hilbcrt gravitational 
Lagrangian density modified by parity- violating Pontryagin density, namely 57c«/3 A ^"'^ ■ It is well-known that it is an 
exact differential, namely, 

Tr[n2] = fi^^ A fl°'^ = d{uja0 A duj°'^ + A A lj-""), (109) 

and therefore it does not yield any field equations as it stands. However, in the model of Jackiw-Pi [65], which is 
based on the Lagrangian density 

= A*r'' + (9f2„/3Af2"'', (110) 

where the term couples to a dynamical scalar field 9, it does contribute to the Einstein field equations. The 

scalar field 9 can be regarded as Lagrange multiplier 0-form imposing the constraint fla^ A Cl"^ = on the resulting 
metric field equations, 

- *G" + 4Di0{de A f2^") -9" A Di^i^{d9 A fl"") = 0, (111) 

which can be found by using (15) and (16). Since the gravitational Chern-Simons term does not contain Hodge dual 
operator, it does not contribute to the coframe field equation directly but through Lagrange multiplier term in (15). 
Compared to Lagrange multiplier terms in the field equations for actions containing quadratic curvature invariants 
studied above, the resulting metric field equations are 3'"'' order partial differential equations for the Chern-Simons 
modified theory. Using the general formula (17) together with the first Bianchi identity and + ^i3a = 0, it is 
possible to show that neither of the Lagrange multiplier terms on the right hand side of (111) has non-zero trace and 
therefore the vacuum field equations require R — 0. 

First order theory based on the Chern-Simons modified gravity (110), with independent coframe and connection 
1-forms has also been studied recently and leads to dynamical torsion [68] . 

In the context of the generalization of gravitational actions we consider, instead of coupling to a non geometrical 
so-called cosmological scalar field 9 it is possible to introduce a term of type f{T) * 1 where *T = fiai3 Ai}"^ is defined 
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for convenience. Therefore, within the same framework used for the other modified models above, it is natural to 
consider the Chern-Simons modified Lagrangian density of the form 

4^jp = 5"a/3 A + /(T) * 1, (112) 

where / = /(T) is an algebraic function of the topological invariant T- 

In order to find the metric field equations for Cgjp, we introduce as for £jp, the total Lagrangian C^^^ = Cgjp + 
I^LM- Using (12), the total variational derivative of the total Lagrangian density can be found as 

^6^^^ = 5ujc.p A {2D{f'n"^) - i(r A - 6"^ A A")} + <5A„ A 6" 

+ ,5r A [inf^" A + u'T - /) *e^ + D\^]. (113) 

Solving for the Lagrange multiplier 2-forms in terms of the other fields and then inserting them into the equations for 
coframe 1-forms we obtain the metric field equations as 

* G" = (/'r - /) * 6*" + 4Dip{df' A il^^") - 6*" A Di^i^{df A O"''), (114) 

with the condition that 0" = 0. In the generalization considered here, the metric field equations (114) remain to be 
3'"'^ order in metric components. Note however that, for the simple choice /(T) = T™, the first term on the right 
hand side is like a cosmological constant proportional to 7"™ and the Lagrange multiplier part, or Z5A" part, contains 
dT instead of the gradient of non-geometrical scalar field 9. In contrast to the equations of motion of the gravitational 
model of Jackiw-Pi (111), the first term on the right hand side appears to be an improvement term. Because of this 
crucial term, the zero trace [R = 0) condition for the vacuum metric equations and the original topological constraint 
*7" = f]"/^ A = for the original metric equations for Cjp are to be replaced with 

i? = 4(/-/'r). (115) 

This equation can be found by tracing the first term on the right hand side in (114), where the Cotton part does 
not contribute to the trace as in the original Chern-Simons gravitational model of Jackiw-Pi. The new "constraint 
equation" (115) allows non-zero values of both of the scalars T and R. The scalar T can be expressed in terms of 
the components of the Weyl tensor relative to the complex null coframe of the standard Newman- Penrose formalism. 
Explicitly in terms of the complex Weyl spinors k with fc = 0, 1, 2, 3, 4, the scalar T can be written out as 

r = 4i {3(*2 _ ^,1) + _ ^-0^-4) + 4(*i*3 - *1*3)} , (116) 

where a bar over a quantity denotes complex conjugation. For Petrov type II spaces, it is possible to find a complex 
null coframe such that = 4*1 = with the other Weyl spinors are non-vanishing. Thus, Petrov type II solutions 
are not allowed in the theory with T = 0. For black hole solutions, which are all of Petrov type D, there is a preferred 
complex null coframe for which the only non- vanishing Weyl spinor is ^2- In this case, T = requires ^'2 to be 
real and therefore, the Petrov D types with real ^2, such as Schwarzschild space-time, are allowed whereas the Kerr 
solution with complex ^'2 is ruled out. Finally, Petrov types III, N and O (conformally flat) one has T = identically. 

VI. THREE DIMENSIONS 

In this last section, we will apply the metric method presented above to the three dimensional cases, in particular to 

the massive gravity models involving modified gravitational Lagrangian densities. The presentation in the section also 
reflects the fact that spacetime dimension n is a crucial parameter in modified gravitational models and demonstrates 
that study of the various modified gravitational actions on a case by case basis can be advantageous. 

In three spacetime dimensions the Einstein-Hilbert action (18) does not lead to a dynamical spacetime model [69]. 
Therefore, various higher order models, such as, third order parity-violating Topologically massive gravity (TMG) 
[70] theory and New massive gravity (NMG) theory proposed recently [71]. At the linearized level about Minkowski 
background, NMG can be considered as an extension of Pauli-Fierz theory to massive spin-2 particle in three spacetime 
dimensions. In a perturbative approach to quantum gravity, NMG has the remarkable properties that it is parity 
invariant, unitary at tree level, ghost- free [71] and renormalizable [72]. Various classical solutions to NMG theory has 
been studied. AdS-wave [73], Kundt wave [74] and black hole solutions [75, 76] to NMG has also been worked out. 
A method, based on Killing vectors, is introduced to study the exact solutions to three dimensional massive gravity 
theories [77] . In this section we study the field equations of NMG theory and express it in a new form in terms of 
Cotton 2-form. We will work with General massive gravity (GMG) theory which is based on the combination of TMG 
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and NMG theories complementing the usual Einstein-Hilbert term [71]. Before we investigate the field equations for 
General Massive Gravity Theory, we first recall some geometrical identities for the curvature tensors which are special 
to three spacetime dimensions. In three dimensions, because the Weyl tensor vanishes identically, the Riemanian 
curvature can be written in terms of contractions and in this case, the identity (88) reduces to the form 

Consequently, the contraction of this identity side by side leads to another useful identity 

A*n"^ = 2R" A*Ra,- ^R'^ *1, (118) 

which will be exploited below. First and the foremost, as a result of the identity (118), in three spacetime dimensions, 
the most general Lagrangian density involving quadratic curvature terms can be written either in the form 



or equivalently in the alternate form 



aR" A*Rc + bR'^ *1, (119) 



aQais A + (26 + |a)i?2 * 1, (120) 



with a, b being arbitrary constants. As we have noted before, this is also the case for the gravitational Lagrangians 

with quadratic curvature terms in four spacetime dimensions as a result of a different geometrical identity (cf. Eqn 
(65) in four dimensions). As a simple cross-check for the various field equations derived in the previous sections, it is 
possible to obtain the metric field equations for the Lagrangian density 

) = if]„^ A + ^ 1^ (121) 

from the metric field equations (68) of the Ricci-squarcd Lagrangian density 

4^) = i?" A *Ra, (122) 

and compare the resulting equations with those that follow directly from (122). In order to do so, inserting the 
appropriate expressions from Eqns. (117) and (118) into the metric field equations *Va = (see Eqn. (68) above), 
one obtains precisely the metric field equations corresponding to the Lagrangian density (121), which is expected as 
a result of the equivalence of the Lagrangians provided by the identity in (118) in three spacetime^ dimensions. These 
considerations are all in accordance with the consistency of the metric field equations for both of the Lagrangians 
(121) and (122) which can readily be written down using the general formulae found in the previous sections, namely 
Eqns. (26), (68) and (81). 

NMG also complements the usual Einstein-Hilbert term with a particular combination of the Ricci-squared and 

scalar curvature-squared terms, which is called X-combination [71]. GMG is an extension of NMG Lagrangian with 
a Lorentz Chern-Simons term and a cosmological constant A, and its original form [71] can be expressed in terms of 
differential forms as 



(3) 

cs 



= A*1-^R*1 + ^(i?« A*i?„ - * 1) + ^ (^"/3 Arfa;^„ + ^uj"^ A A uj\) , (123) 

where /i is the mass parameter in TMG theory whereas rri^ is the mass parameter in NMG theory. Therefore, the 
GMG Lagrangian, which has two mass parameters, reduces to that of TMG Lagrangian in the limit m oo and 
to that of NMG Lagrangian in the limit fi oo. Note that the GMG Lagrangian density (123) can equivalently be 
written in the form 

4mg = a * 1 - ii? * 1 + j^in^p A - ii?2 ^ 1) + ^ (^"^ A duj^^ + A ujf'^ A uj\) , (124) 

by using the identity given in (118). The variational derivative of the Lagrangian density (123) (or equivalently, (124)) 
and the resulting Lagrange multiplier 1-forms have special properties that deserve further scrutiny. We choose to work 
with the form (124) and the constrained total variation derivative of the Lagrangian density Ctot = ^gmg + ^lm 
leads to 

SCtot = StJae A [^D * (O"'^ _ l^J^a^) + _ Ijj ^ Qcfi _ 1 A - A A«)} 

+ ^^aA{(l + ^)*G« + (A+^)*r + 5^*T«[0]+m«}+a„Ae«, 
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(up to an omitted closed form) where the definition of the canonical energy momentum 2-form for curvature 2-forms 
defined in Eqn. (80) has been used for convenience as before. Note here that the trace of the energy momentum 
2-form (80) in three dimensions is 

T«„[0] * 1 = -ifi^^ A (125) 

The Lagrange multiplier 1-form can be solved in terms of other fields using the field equations obtained from the 
connection variation using general formulae (16) and (117). This yields 

= :^iaD * (0"'5 - \R6''^) + ]j^{R^ - \Re^). (126) 

Here, the first term is the contribution of Ck terms and are fourth order in metric components, whereas the last 
term contributes to the metric field equations as the Cotton 2-form C" = D{R" — jR9") and is third order in metric 
components. Although Cotton 2-form is defined in arbitrary dimensions and Cq, (8) 0" is invariant under conformal 
transformation, its vanishing is necessary and sufficient condition for conformal flatness in three dimensions. In 
addition, the Cotton 2-form is traceless, iaC" = 0, satisfies Bianchi type identity, 6a A C" = and C" is covariantly 
constant, DC" = 0, in three dimensions [78]. 

The significance of the particular choice of the Lagrangian density Ck in (124) shows up in the expression for the 
Lagrange multiplier terms, namely DX^ . The first of the terms in (126) is NMG part whereas the second term is the 
TMG part. The traces of each parts vanishes separately, and therefore they do not contribute to the trace of the 
metric field equations of the GMG theory, to put it mathematically one has Oj^ A DX^ = for the Lagrange multiplier 
given in (126). Consequently, as can be explicitly verified from the resulting field equations, the trace of the terms in 
the field equation coming from the Lagrangian density Ck add up to Ck up to a sign. (We refer to the trace formulae 
given in Eqns. (25), (69) and (82)). Using the result (126), the GMG metric field equations can be written as 

(1 + ^R)*^^ + (A + ^i?2) *r + iC« + ^ *T«[f2] 

-4^D*{dRAe") + ^Di0D*n^" = O. (127) 

As we shall exemplify below, the metric field equations given in (127) for the GMG field equations seem to be more 
advantageous compared to the original form of the field equations, since for example, for the fact that it can be also 
expressible in parts or as a whole in terms of Ricci 1-forms and scalar curvatures using (117). Note that the TMG 
theory is a third order theory whereas NMG is a fourth order theory. Using (117) and the fact that in three dimensions 
the Ricci 1-forms and scalar curvature satisfy the identity iaDR" + ^dR = (This identity can be obtained using the 
property that Cotton 2-form is traceless) one finds that the contribution of the Lagrange multiplier term, i.e. DX^ 
term, to the metric field equations takes the form 

£)A" = ;^Z)*C" + iC". (128) 

This form of the Lagrangian multiplier also makes it easier to see the various properties of the NMG theory in 
comparison to TMG. The first term on the right hand side of (128), which is fourth order in metric components, is 
unique with the properties: (i) Trace-free, (ii) Covariantly constant, (iii) Parity-preserving. In contrast, the TMG 
part of the Lagrange multiplier shares these properties except that it is parity-violating. 

It is shown that TMG is "square root" of the NMG theory at the linearized level [71] . The presentation of the field 
equations above makes it possible to sharpen this result. In fact, such a relation between the two theory is encoded 
in the Lagrange multiplier structure of the GMG theory and can be put in mathematically precise form under certain 
assumptions using Eqns. (127) and (128). If we assume for simplicity that A = and i? = 0, then using (127) the 
TMG and the NMG field equations can be written as 

(*L> + Ai)i?" = 0, (129) 
{*D*D + Tn^)R" = -^T"[n], (130) 

respectively. Eqn. (129) can be derived by simply applying Hodge dual operator * to the TMG equations obtained 
from (127). The energy momentum 1-form T"[f2] in the right hand side of Eqn. (130) can be expressed in tems of 
Ricci-squared terms for R = 0. This correspondence amounts to the fact that both TMG and NMG can be written 
in Dirac type and Klein-Gordon type equations respectively. This point of view is elaborated in [79] relative to a 
coordinate coframe where the correspondence is utilized to provide a method to generate exact solutions for the NMG 
theory for any given Petrov-Segre type D and N solutions of the TMG theory [69, 78, 80, 81]. In particular, the 
non-dispersive, line-fronted gravitational wave metric of the Kerr-Schild form 



g = 2H{u, x)du ® du + dv ® du + dv ® du + dx ® dx 



(131) 
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expressed in terms of local coordinates {u, v, x} where u, v are real null coordinates and x is a spatial coordinate. 
The metric (131) linearizes the TMG equation (129) for the real metric function H(u,x) [82]. Moreover, for the 
Pctrov-Segre type N metric (131), one has r"[f7] = identically and therefore, it also linearizes the NMG equations 
(130). Consequently, there is a class of gravitational wave solutions of the form (131) which is common to both TMG 
and NMG theories. 



VII. CONCLUDING REMARKS 



We have presented a unified and practical method to obtain field equations for some modified gravitational La- 
grangians containing quadratic curvature invariants and their algebraic generalizations in n ^ 3 spacetime dimensions. 
We have adopted the language of differential forms throughout our presentation which provides further insights even 
for the well-known results. As clearly indicated by our presentation for the GMG theory in three dimensions relative 
to an orthonormal coframe, the derivation of the metric field equations and the Lagrange multiplier terms provides 
alternative forms for the field equations and insight for the striicture of them. 

Although we did not address the question of coupling matter to specific gravitational model, our treatment can be 
even extended to the case to in which spinor matter fields present or to the case where non-metricity is present. We 
presented the variational derivatives with respect to independent coframe and connection 1-forms in an explicit form 
for all the generalized gravitational actions such that these extensions can easily be pursued. The presentation also 
makes it possible to compare the Palatini formalism with the metric formalism. Rather then providing an exhaustive 
list of modified gravitational Lagrangians and their corresponding field equations we presented a unified method 
of the derivation for the field equations using popular gravitational models. We leave the further studies such as 
the exact solutions, and other physical and mathematical features of the new models we put forward, i.e., modified 
Chern-Simons theory and GMG theory to forthcoming research. 



Appendix 

In this section we will briefly explain the notation. For further details we refer to, for example, [27] and [47]. The 
definitions of basic geometrical quantities we use can be summarized as follows. The metric relative to an orthonormal 
coframe 1-forms 0" is g = rjafiO" 0^ , the wedge products of basis 1-forms arc abbreviated as 0°"^'" = 9" A 9^ A ■ • ■ . 
The contraction with respect to frame fields is denoted by where iad'^ = Sp. Oriented volume element can be 
written as *1 = 0Oi - (n-i) ■y^,JJg]-g ^ jg linear Hodge dual operator and ** = (— 1)''+p("^p) acting on a p-form in an 
n dimensional spacetime (s is the metric signature). The indices of the permutation symbol Saia2---ctn relative to an 
orthonormal basis is raised and lowered by Tjais = rj"^ and £oi23---(n-i) = +1- Cartan's first structure equations can 
be used to define torsion 2-form and torsion tensor as 



1-fo 

2-forms in terms of Riemann tensor are 



where a;"^ = oj'^^9'^ are connection 1-forms. The contorsion 1-form is defined as 6" = K"^ A 9^. The curvature 



(Cartan's second structure equations). Ricci 1-forms can be defined as the contraction of the curvature 2-forms 



where R^x are the components of the Ricci tensor. The scalar curvarture is i? = iaR". Similarly, the components of 
the Weyl 2-form and those of the Weyl tensor C"^^^ are related by C"^ = ^C^p^J^". 

A tensor cj) of type (r -|- p, s) which is totally antisymmetric with respect to p number of indices can uniquely 
be associated with the (r, s) tensor-valued p-form 4'°'^"^ p-^^p^-'-Pr' -R-elative to an arbitrary coframe considered as 
multilinear maps on (co)framc fields they arc related by 

^"'"'■■■"''/3ife.../3,.('5/3i, • • • , e^J = 0(0"! , 6'"^ ■ • ■ , 61"" , e^, , e^^, • • • , e^j^ , e^, , e^j^ , • • • e^J. 

Covariant exterior derivative D acts on tensor- valued p- forms to yield tensor- valued (p -|- l)-form preserving the 
tensorial type as follows 

/3i/32-/3r - /31/32-/3. A A0 /3i/32-/3. /3i A A/32-/3. 
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More generally, for tensor valued p-form (j)^ and tensor valued g-form ip"^ 

Acting on scalars D reduces to exterior derivative d. The expression _D(/)"^"^ '""^^^^...^ can easily be related the 
covariant derivative V x4> using the relation between the exterior and Riemannian covariant derivatives, namely d = 
6'^ A Va- The covariant derivative = X"Va can be defined by its action on basis coframe: Vx^" = —{ixi^^s)^^- 
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